Sudden Collapse of a Granular Cluster 
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Single clusters in a vibro-fluidized granular gas in connected compartments become unstable 
at strong shaking. They are experimentally shown to collapse very abruptly. The observed cluster 
lifetime (as a function of the driving intensity) is analytically calculated within a flux model, making 
use of the self-similarity of the process. After collapse, the cluster diffuses out into the uniform 
distribution in a self-similar way, with an anomalous diffusion exponent 1/3. 

PACS numbers: 45.70.-n, 05.45.-a, 02.60.Lj 



One of the key features of a granular gas, making it 
fundamentally different from ordinary molecular gases, 
is its tendency to spontaneously separate into dense and 
dilute regions Q. This clustering originates from the 
dissipative nature of the particle collisions. It is an un- 
wanted effect in many applications where granular mate- 
rial is brought into motion. Therefore we study (within a 
simple geometry) how declustering occurs. We find that 
the breakdown of a cluster can be very abrupt, making 
declustering very different from clustering in reverse time 
order. 

The experimental system (see Fig. 1) consists of a row 
of N equal compartments, separated by walls of height 
h = 25.0 mm and filled with a few hundred steel beads of 
diameter 3.0 mm. We start out with all the particles in 
the middle compartment and bring them into a gaseous 
state by shaking the system vertically. For weak shaking 
the cluster is stable: after some initial spilling, a dynami- 
cal equilibrium is established between the outflux of slow 
particles from the cluster and the influx of fast particles 
from outside ||, ^, ||. For sufficiently strong shaking, 
however, the cluster breaks down. 

Two different regimes are observed: (i) At very strong 
shaking the breakdown occurs immediately, and the clus- 
ter spreads out over the boxes with its profile widening 
as t^^"^ (instead of the standard t^/^ diffusion law), (ii) 
At less violent shaking, the cluster seems to remain sta- 
ble for a long time, showing only a small leakage to its 
neighbors. But suddenly it collapses and subsequently 
diffuses over all boxes. The sudden death of the cluster 
is depicted in Fig. |l|. 

One thing this figure shows is that the breakdown of a 
cluster is strikingly different from the reverse process of 
cluster formation, which is known to take place gradually 
and (for all > 3) via transient states showing clusters 
in several boxes ^. This lack of time-reversibility is 
yet another consequence of the dissipation in the system. 

The abruptness of the collapse allows us to define a 
cluster lifetime r (via nci{T) = 0, see Fig. 0e). In Fig. | 
the measured lifetimes are plotted as a function of the 
inverse shaking strength B (Eq. ||) for various values 
of TV. The data lie on a universal envelope curve, until 
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FIG. 1: Four images from a 5-box experiment, at driving 
parameter B — 8.2. The cluster is clearly present until t = 42 
s, then suddenly collapses, leaving no trace one second later. 
Plot (e) shows the time evolution of the cluster fraction nci{t), 
evaluated from the flux model. 
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All the above experimental observations can be ac- 
counted for quantitatively by means of the flux model of 
refs. 1^, 1^, ^. At the heart of this model is a flux function 
F(nfc), describing the outflow from the fcth box to each of 
its neighbors. It is a non-monotonic function of nj. (the 
particle fraction in the box): F{nk) first increases with 
rife, but beyond a certain value of Uk it decreases again, 
as the increasingly frequent inelastic collisions slow the 
particles down, so that they cannot make it over the wall 
to the neighboring compartments anymore. The precise 
form of F{nk) is not very critical, as long as it is a one- 
humped function. We will use Eggers' flux function |^ 



with 



F{nk) = CVBnle'^" 
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FIG. 2: (a) Lifetime r vs driving parameter B, determined 
from experiments with A*' = 3, 5 boxes (dots with error 
bars), and from numerical evaluation of the flux model with 
N = 3,5,10,20,40,80 boxes (empty symbols). The solid 
curves are analytical solutions for A'^ = 3 Q and for the enve- 
lope curve, which goes roughly as exp((3/2)B) (cf. Eq. 
(b) Enlargement of (a), showing the experimental results for 
A'' = 3 and A'' = 5 in more detail. Every point is based on 
15 repetitions of the experiment; the vertical error bars de- 
note the maximal deviation from the average r measured, the 
horizontal ones represent the accuracy in B. (c) Bifurcation 
diagram for N — 5, showing the critical value Be, 5. Declus- 
tering occurs to the left of this value. 



The driving parameter B depends on the total number 
of particles P and their properties (radius r, restitution 
coefficient e of the particle collisions) , on the geometry of 
the system (height h of the walls, ground area S of each 
box), and on the frequency / and amplitude a of the 
shaking. The factor C only determines the absolute rate 
of the flux, and can be incorporated in the time scale. 

The equation of motion for the fraction in box A: is 0] 



2F{nk)+F{nk+i), 



(3) 



where k = 1,2,..,A^. Here we assume a nearest neigh- 
bor interaction, and a cyclic arrangement of the boxes 
[k = N + 1 equals fc = 1). We further impose particle 
conservation, X]fc"-fc = 1- 

The numerical results shown in Fig. |l]e and Fig. || have 
been obtained using the above flux model, starting out 
with all particles in one box (labeled cl). They quanti- 
tatively agree with the experimental observations. The 
decaying cluster goes through three different stages: 



The starting stage is a very short one, in which both 
Jici and F{nci) display a jump compared to Ui and F{ni) 
in the surrounding boxes, i — 1,2, .. (we have to consider 
one side only because of the symmetry in the system). 

In the second stage, the flux has become continuous 
but the particle fraction remains discontinuous. How- 
ever, its low-density counterpart no (defined by F{no) = 
Firici)) does continuously connect to ni. We will use 
this fact later in the analysis of the envelope curve. The 
flux gradually grows, and eventually F{nci) reaches its 
maximum value. This is accompanied by rapid density 
changes and the sudden death of the cluster at the life- 
time r. 

In the third and last stage, both rici and F{nci) flt 
continuously to the other boxes (see Fig. ^). The 
remains of the cluster diffuse over the whole system until 
the uniform distribution is reached. 

In what follows we will analytically solve the flux 
model. First we focus on the third stage. We rewrite the 
problem into its continuum version, by setting n{x, t) = 
n-kit) {x = kw by definition, where the box width w will 
be incorporated in the x-scale). Eq. ^then becomes: 



dtn = da:xF{n{x,t)) 
^C\/Bd,, (n{x,t) 

condition 



(4) 



takes the form 



and the conservation 
J^^n{x,t)dx = 1. 

For very strong shaking (regime (i), where t is van- 
ishingly small) the diffusive stage sets in almost immedi- 
ately. Here i? — > 0, and Eq. H reduces to 



dtn = CVBdxxin') = 2CV B{{dxny + n9,,n), (5) 

which is known as the porous media equation ^, |^ . The 
decay of the cluster in this limit is depicted in Fig. It 
is self-similar: all curves in Fig. ^ fall onto a single curve 
if we properly rescale the axes (Fig. ^b). The original 
partial differential equation (PDE) can thus be brought 
back to an ordinary differential equation (ODE) in terms 
of the self-similarity variable rj — x/{CB^/^ty^'^. With 
n(x,t) = H{ri)/{CB^/'^tY/^, Eq. | now takes the form: 



0. 



(6) 



Its symmetric solution is H{rj) ~ Hq — (1/12)?7^ (with 
the constant Hq ~ (3^/'^)/4 « 0.361 determined by 
H{ri)dr] = 1). This inverted parabola, depicted in 
Fig. ||b, represents in one curve all the stages of Fig. ||a. 
The scaling of the axes shows that the height of the clus- 
ter decreases as t~^^^, and its width grows as t^^'^. This 
anomalous diffusion (with exponent 1/3) is also found 
in porous media |^]. The slowed down diffusion of 
the front originates from the quadratic n-dependence in 
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FIG. 3: (a) The diffusing profile at successive times t (in 
the limit B = Q). (b) The function H{rii) onto which all the 
curves in plot (a) collapse by a proper rescaling of the axes 
{H oc nt^^^ and oc xt~'^'''^). 



dxxn^ (Eq. ||). In the context of our granular model this 
reflects that particles only diffuse to neighboring boxes 
through 2-particle collisions: In the strong shaking limit 
the absence of particles slows dovifn further diffusion, the 
presence enhances it. 

For non-zero B (regime (ii), where the diffusive stage 
has to wait until after the sudden death) we have an 
additional dimensionless variable x = B/C{t — r). 
Its influence diminishes with time and the solutions 
converge to the inverted parabola of the case _B ^ 0. 

Next we turn to the second, semi-continuous stage. For 
moderate shaking, B « Bc^n, this stage can take quite 
a long time. At the critical point Bc,n the lifetime r 
even diverges to infinity and the cluster becomes stable. 
Starting out from the initial ..,0,0,1,0,0,.. distribution, 
the system first very slowly approaches a distribution in 
which Uci is close to the cluster density at the saddle-node 
bifurcation (see Fig. ^) and all other boxes contain equal 
fractions n/, = {l — rici) / {N — 1). Only when it has passed 
this phantom equilibrium (i.e., when rid is below the level 
of the saddle-node bifurcation) , the system quickens its 
pace and the sudden collapse occurs. 

This means that t is the time it takes to pass the 
phantom equilibrium. It can be calculated either 
numerically or analytically (by integrating the Taylor 
expansion of Eq. ^), with the result r oc (-Bc.w ^ B)^^/^. 
So T diverges as the inverse square root of the distance 
to the critical point, which is the common (mean field) 
power-law behavior near a second order phase transition 
as we have here for r 

Finally, we calculate the lifetime away from the crit- 
ical point, i.e., for B considerably smaller than B^^n- 
This will give us an analytical expression for the enve- 
lope curve in Fig. For these B values, the collapse 
occurs before the particles leaking out of the cluster have 
had time to fill the outermost boxes to any significant 
level. Therefore, the behavior does not depend on the 
value of N: The system does not feel its finite size dur- 
ing the cluster's lifetime, so the number of boxes can be 



taken to be infinite. 

The time-evolution of the cluster is described by Eq. |^, 
with F{n_i) = F{ni): 



duci 
~dt 



= -2F{nd) + 2F{ni) = ~2F{no) + 2F{ni). (7) 



This equation contains rii, which is governed by a simi- 
lar equation of motion (Eq. ^ containing n2, etc. So we 
have to deal with an infinite number of coupled nonlin- 
ear ODE's This is a problem that cannot be solved 
directly, so we attack it in five steps. 

Step 1 : We first rewrite the problem into its continuum 
version, and replace the cluster density rid by its low- 
density counterpart n(0,i) = no(i) ~ nc/exp(— _Bn^j/2). 
Thus, without influencing the fluxes (since F{nQ) — 
F{nd)) we make n{x, t) continuous in x = 0. The density 
n{x, t) obeys Eq. ^ plus a conservation condition saying 
that the increase of material into the rest of the system 
equals the influx from x — 0: 



dt n{x,t)dx^-{dx{F{n{x,t))) (8) 
Step 2: The cluster-equation Eq. |^ now becomes: 



dn. 



2(dxF{n{x,t))"^ 
= 2F'{n{Q,t)){dM^,t))x=o- 



x=0 



(9) 



Since F' can be derived directly from Eq. the problem 
reduces to determining dxn{x,t) at a; = 0. 

Step 3: In order to do so, we observe that changes in 
n(0,i) happen on a much longer timescale than in the 
surrounding boxes , so the cluster acts as a constant 
reservoir spilling granular material. This approximation 
is illustrated in Fig. ^: the profile in the system builds 
up while no remains constant. In fact, this build-up takes 
place in a self-similar way (see Fig. |b). So, as before, 
the problem for n{x,t) can be formulated in terms of one 
variable ^ = x/{noCB'^/HY/^ . Setting n{x,t) = noG{C), 
Eq. ^ becomes an ODE for G'(^), and also the accompa- 
nying conservation condition (Eq. |^) contains ^ only: 



p/3 



«=o 4(1 - f3) 7o 



GiOd^, 



(10) 



Bnl 



where /3 = 

Step 4'- The slope {dxn{x,t))x=a can now be approx- 
imated by riQ(t){d^G{S,))^=Q{dxS,)x=o, where we have re- 
vived the (slow) time-dependence in no{t). With dxS, = 
{noit)CB^/H)-^/^, Eq. I becomes: 



~~dt 



= 2^^F'{n,{t)){d,G{0),^,. 



(11) 
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FIG. 4: (a) The density profile in the boxes surrounding 
the cluster, at 20 consecutive (equidistant) moments in time, 
within the constant no approximation; no = 0.223 and B — 
3.00. (b) The function G = n/no onto which all curves in plot 
(a) collapse by a proper rescaling of the axes oc xt~^^^.) 



All quantities on the right hand side are tractable. We 
re-express no{t) in terms of nci{t), derive F' from Eq. 0, 
and solve Eq. ^ to determine (9^G(^))^^. This last 
step still requires some work, because Eq. |lO| does not al- 
low an analytical solution and moreover contains (3 (and 
hence no(i)) explicitly. Since /3 is small, however, we may 
expand G(^) and Eq. |l^ in terms of /S, and solve numer- 
ically. In leading order we find [d^G{^)) = —K = 
—0.3138. Inserting all this, Eq. pT| becomes: 



(12) 



This is an ODE for rici in closed form, which replaces the 
original problem (Eq. |^) consisting of an infinite number 
of coupled ODE's. 

Step 5: Finally, we integrate Eq. |l^ over the cluster 
density (dropping the suffix cl) and find an analytical 
expression for the lifetime r away from the critical points: 



1 v/cBi/2^-3/2e(3/4)^"'rfri 
8i^(l - Bn2e-^"') 



(13) 



Here rithr is the value of n at which the sudden death 
occurs. For the evaluation of the solid curve in Fig. ^a 
we used nthr — 0.5, but this value is not too critical 
(cf. Fig. |]e). The only free parameter is the constant C: 
if this is adjusted properly, the analytical r curve agrees 
with the measured data over the whole range of B-values. 

The above expression shows that t roughly increases 
as exp(i?'^/2). Recalling that B is the inverse shaking 
strength, this underlines the experimental observation 
that even a small reduction in the shaking strength 
causes a tremendous increase of the cluster lifetime. 



In conclusion, in the studied compartmentalized 
system clusters break down very abruptly, in contrast 
to their slow formation. As clustering itself, the lack 
of time reversibility originates from the dissipative 
nature of the particle collisions: The breakdown of the 
unstable cluster is delayed because most of the energy 
input is dissipated through collisions in the cluster. The 
dynamics is quantitatively described by a remarkably 
simple flux model, which can be analytically solved. 
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The case A'^ = 3 is contained in a single ODE, which can 
be solved exactly: r = dn/2{F{(l - n)/2) - F{n)} 
(solid curve in Fig. 2), with nthr determined from n(r) = 
0. An analogous expression holds for general A'^ close to 
the critical point, where all the surrounding boxes are 
flUed to the same level Uk ~ {1 — nci)/{N — 1). 
A cluster change And (flowing into the surrounding 
boxes) corresponds to a much smaller change Ano 
in its low-density counterpart: dud/dno = (1 — 
Bn2,)-ie^"c,/2 « (1 _ _B)-^e^/2 ^-^^^ ^ ^ 

most B-values of interest (excluding the strong-shaking 
regime B ^ 0) \Anci \ = J2i^o ^ ^'^o- 



